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Abstract 

Even though the KdV and modified KdV equations are nonlinear, we show that suitable linear combi¬ 
nations of known periodic solutions involving Jacobi elliptic functions yield a large class of additional 
solutions. This procedure works by virtue of some remarkable new identities satisfied by the elliptic 
functions. 
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The Korteweg-de Vries (KdV) and modified Korteweg-de Vries (mKdV) equations are known to 
possess many remarkable properties. They are examples of completely integrable systems with soliton 
solutions which have diverse applications of physical interest. Soliton solutions of the KdV equation 

VjI QuUx T Uxxx — 0 ; (f) 

are discussed in many texts ||l|. The simplest, periodic, cnoidal travelling wave solution is 

ui{x, t) = — 2Q!^dn^(^i, m) + /Sa^ , = a{x — bia'^t) , (2) 


where a, m and (5 are constants, and the “velocity” bi is given by 


6i = 8 — 4m — 6/3 . 


(3) 


Here, we use the standard notation dn(^,m), sn(^,m), cn(^,m) to denote Jacobi elliptic functions, 
where m is the elliptic modulus parameter (0 < m < 1). Note that the period of solution (2) 
is 2K{m)/a, where K(rn) is the complete elliptic integral of the first kind In the limiting case 
m = 1, one recovers the familiar result —2Q:^sech^(a(x—6ia^t)). The properties of periodic multisoliton 
solutions more general than eq. (2) have also been discussed ||3|. 

In this paper, we make suitable linear combinations of solution (2) to obtain a large class of new 
periodic solutions of the nonlinear KdV equation. This procedure works by virtue of some remarkable 
identities satisfied by the Jacobi elliptic functions, which to the best of our knowledge, do not seem to 
be previously known. Our solutions consist of adding terms of the kind given in (2) but centered at p 
equally spaced points along the period 2K[m) ja, where p is any integer. The p-point solution is 

Up{x, t) = —2c^ ^ + (io? ; di = dn[^p H— — -^m] , = a{x — bpO^t) . (4) 

i=i P 

Clearly, p = 1 is the original solution, but for any other p, we have new solutions of period 2K{m)/pa. 
For convenience, we also define the quantities Si and c* in analogy to the quantity di defined above: 






m] 


2{i-l)K{m) 

a = cn[4p H-, ml 


p 


(5) 


To see why eq. (Q) is a solution of the KdV equation (||), it is best to first consider several small 
values of p. For p = 2, substitution of eq. (Q) into the left hand side of the KdV equation (|^ gives on 
simplification: 


4ma^{(8 - 4m - 62 - 6/3)(siCidi + 5202 ^ 2 ) + l2{d\s2C2d2 + d\sicidi)]. ( 6 ) 
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One gets a solution only if the above expression vanishes. For p = 2, the quantities di = dn(,f 2 )W') 
and d 2 = dn[^2 + K(m),m] are well-known to satisfy = 1 — m, which on differentiation 

gives df 5202^2 -|- d^sicidi = 0. Therefore, U 2 {x,t) is a solution of the KdV equation with velocity 
62 = 8 — Am — 6P, the same velocity as 61 . 

For p = 3, substitution of eq. (|^ into the left hand side of the KdV equation (^) leads to the 
expression: 


4ma^{(8 - 4m - 63 - 6/3){siCidi + S 2 C 2 d 2 + 5303 ^ 3 ) 


12[d‘l{s2C2d2 + sscsds) + dl{s^c^d^ -|- sicidi) d|(sicidi 5202 ^ 2 )]}- 


(7) 


In this case, di = dn(^ 3 ,m),d 2 = dn (^3 -|- 2iF(m)/3, m), ^3 = dn (^3 -|- AK{m)/3,m). In contrast to 
the p = 2 case, there are no well-known relations between Jacobi elliptic functions at the arguments 
chosen above. However, we have recently discovered numerous new identities @ which reduce by 2 
(or a larger even number) the degree of cyclic homogeneous polynomials in Jacobi elliptic functions. 
The relevant identity for the p = 3 case is 

d‘l{s2C2d2 + S3C3d3)+dl{s3C3d3 + SiCidi)+dl{siCidi + S2C2d2) = ('SlCldl-h 5202^2+ S3C3d3) , 

( 8 ) 

where 

Q = sii^[2K{m)/3,m] . (9) 

Identity (^) can be established analytically using the addition theorems @ for Jacobi elliptic functions. 
It is also easy to check numerically. Substituting identity (^) into expression ( 0 ). shows that U 3 {x,t) 
is a solution of the KdV equation with velocity 

( 10 ) 


( 11 ) 


63 = 8 - 4m - 6/3 -k 24 ( 1 - - 


Similarly, the identity analogous to (0) for p = 4 is 


44 4 

df SjCjdj = — 2 \/l — m SiCidi 
i=l j^i 


2 = 1 


and consequently, the expression for the velocity is 

64 = 8 - 4m - 6/3 - 2AVT — m . 


( 12 ) 


For an arbitrary integer p, we have found that the relevant identity is 

p p p 

2 = 1 j^i 2=1 


(13) 
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where the constant A{p, m) can be evaluated in general by choosing any specific convenient value of 
the argument ^ of the Jacobi elliptic functions. The values at m = 0,1 are particularly simple: 

A{p, m = 0) = —^{p — l){p — 2) ; A{p, m = 1) = 0 . (14) 

The general expression for the velocity of our p-point solution is 

bp = 8 — 4m — 6/3 + 12^(p, m) . (15) 

In Figure 1, we plot the velocity hp as a function of the modulus parameter m for various values of 
p. In making the plots, we have chosen /3 = 0. Note that for this choice of /3, the velocity is always 
positive and corresponds to right moving waves for p <3, whereas the velocity changes sign for p > 3. 

Two types of mKdV equations are discussed in the literature corresponding to the sign chosen 
for the nonlinear term. We will consider these two types of mKdV equations separately, in order to 
stay with real solutions. The first type of mKdV equation is 

Vt - + Vxxx = 0 . (16) 


The best known periodic travelling wave solution is 

vi{x, t) = ±\/ma sn(? 7 i, m) , rji = a{x — qia^t) , 


(17) 


with “velocity” qi = —(1 + m) and period AK{m)/a. In the limiting case m = 1, ui(x,t) reduces to 
the familiar form ±atanh(a(x + 2a^t)). 

For any odd integer p, we find the following solutions of the mKdV equation ([^) by a specific 
linear superposition of the basic solution (|T^: 

p 

Vp{x, t) = ±\/rna Sj , p = odd , (18) 

i=l 


where we define 


Si 

di 


sn[r]p + 
dn[r]p + 


4(i — l)K(m) 
P 

4(z - l)K{m) 

p 


m] ; Ci = cn[r]p + 


m] ; rip = a{x — 


4(z - l)K(m) 

p 

qpo'^t) . 


,m] ; 


In order to verify that Vp{x,t) satisfies eq. (^), one needs the identity 


p p ^ 

m ^ ^ Cjdj + 2m 

1 ___ 

' P 

^ ^ ^kdk 

= {B{p, m) — C{p, m)} 

' p 

^ ^ ^kdk 

i=l j^i 

_i<j 

.k=l 


.k=l 


(19) 


( 20 ) 
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which can be established from the much simpler identities y] 

p p p 

m'^SiSj = B{p,m) ; m ^ SiSjSk = C{p,m)'^Si , p = odd, 
i<j i<j<k 2=1 

by using the properties of Jacobi elliptic functions. The general expression for the velocity is 

Qp = -{1 + m) - 6{B{p,m) - C{p,m)} . 

Some explicitly computed values of the constants B{p,m) and C{p,m) are: 

B{l,m) = 0 ; C{l,m) = 0 ; B{3,m) = —mQ ; C{3,m) = —IjQ . 


(21) 


( 22 ) 


For any even integer p, the linear superposition of elementary solutions does not work, since identity 
( ^iD is restricted to odd values of p. However, remarkably enough, for this case, we find that products 
of elementary solutions do give new solutions! For example, the solutions for p = 2 and p = A are: 

V2{x,t) = AzamsiS2 ; Vi{x,t) = ±Q;m(l — y/l — m)siS2SsS4 . (23) 


These solutions can be verified using the general identity Q 

CidiCjdjCkdk , s Cidi 

y -= L>(p, m) y - , p = even 

. , SiSiSk ^ Si 

i<j<k ® *=i * 

The expressions for the velocities are 


(24) 


q 2 = —2(2 — m) ; qa = —2(2 — m) — 12VT — m . 


(25) 


It is well-known [Ql that if n(x, t) is a solution of the mKdV equation (pl^), then the Miura transform 
u{x,t) = ± fa; is a solution of the KdV equation (||). Therefore, using solutions (^) and (^), we 

can obtain even more solutions of the KdV equation via the Miura transform. The simplest solution 
{p = 1) coming from eq. (O) is 


u{x,t) = a‘^[msn‘^{ri,m) ± \/mcn(r/, m)dn(? 7 , m)] 


(26) 


with velocity —(1 -|-m), corresponding to a left moving wave. In view of the identities (^), it is clear 
that for odd integer p, appropriate linear superposition of the simplest solution (|^) also solves the 
KdV equation with the same velocity qp given in eq. (|^). Similarly, the solution resulting from a 
Miura transform of V 2 {x,t) in eq. (^) is 

u{x,t) = 0 ^ 1171^8182 Az m{siC 2 d 2 + S 2 Cidi)] . (27) 
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However, this is not a new solution since it simplifies to give either Q;^(2msf — m) or a^{2ms2 — m), 
both of which are KdV solutions (^) with (3 = 2 — m. 

Finally, we consider the second type of mKdV equation 

Vt + Qv^Vx + Vxxx = 0 . (28) 

Here again, we have found several new solutions via linear superposition: 

p 

Vp{x, t) = ±a'^di , qp = 2-m + Q{E{p, m) - F{p, m)) , {p = l, 2, 3,...), (29) 

i=l 

P 

Vp{x, t) = Ey/ma ^ Ci , Qp = 2m — 1 + 6{G{p, m) — H{p,m)) , p = odd, (30) 

i=l 

p-i 

Vp{x,t) = ±a ^ {di - di+i) ,qp = 2 - m- 6{I{p,m) - J{p,m) + L{p,m)) , p = even, (31) 

i odd 

where the constants E{p,m), F{p,m), etc. are defined by the identities [Q] 
p p p 

didj = E{p, m) ; Y didjdk = F{p, m)Ydi ] 
i<j i<j<k i 

P P P 

mY^iCj = G{p,m) ; m Y CiCjCk = H{p,m)Y^i 5 
i<j i<j<k i 

P P 

Y didj = I{p, m) , {i+ j) odd ; Y didj = J{p, m) , (i + j) even ; 

i<j i<j 

P P P~^ 

^ ^ didjd}^ ^ ^ didjdj^ — L^P^tti') ^ ^ {di d^+i) • 

i<j<k i<j<k i odd 

i^j+k=odd 2 +j'+fc=even 

Some explicit simple values at m = 0 are: 

E{p,0) =p{p-l)/2 , F{p,0) = {p-l){p-2)/6 , G(p,0) = 0 , H{p,0) = {p^ - l)/6 , 
l{p,0) =p‘^/4 , J{p,0) =p{p-2)/4 , L{p,0) = {p-l){p-2)/6 . (33) 

At m = 1 all constants E, F, G, H, /, J, L vanish for any p. Also, the constants for arbitrary m but 
specified small values of p are: 

E{2,m) = I{2,m) = ^ F{2,m) = J{2,m) = L(2,m) = 0 , 

E{3, m) = 1 — mQ + 2-\/l — mQ , F(3, m) = (1 — Q)/Q , 

G{3,m) = —m(l — m)Q/{l — mQ) , H{3,m) = (1 — mQ)IQ , 

£'(4, m) = 2r(l + r + P) , £(4, m) = L(4, m) = P , 

/(4, m) = 2r(l + r^) , J(4, m) = 2p , r = (1 — m)^/^ . (34) 
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Our method of judicious linear superposition for finding new solutions also works for many other 
nonlinear equations. Indeed, we have verified that it works for the nonlinear Schrodinger and KP 
equations as well as the model. These results will be reported in a forthcoming publication |^. 

The solutions we have obtained for the KdV equation all correspond to one gap periodic potentials. 
This process can be generalized to obtain periodic potentials with a finite number of band gaps ii- 
One of us (A.K) thanks the Department of Physics at the University of Illinois at Chicago for 
hospitality. We acknowledge grant support from the U.S. Department of Energy. 
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Figure Caption 

Figure 1: Velocity bp of the periodic KdV travelling wave solution Up{x,t) [eq. (^] as a function of 
the modulus parameter m for p = 1,2, 3,4. 
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